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We consider in the framework of the fluid/gravity correspondence the dynamics of hypersurfaces 
located in the holographic radial direction at r = tq. We prove that these hypersurfaces evolve, 
to all orders in the derivative expansion and including all higher curvature corrections, according 
to the same hydrodynamics equations with identical transport coefflcients. The analysis is carried 
out for normal fluids as well as for superfluids. Consequently, this proves the exactness of the bulk 
viscosity formula derived in larXiv: 1103. 1657. via the null horizon dynamics. 
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The AdS/CFT correspondence [i|, y] relates certain 
strongly coupled large N gauge theories to weakly cou- 
pled classical gravity, thus enabling the calculation of 
strongly coupled field theory quantities via General Rel- 
ativity (GR). In particular, the map of thermal states 
of field theories to classical gravity solutions contain- 
ing black holes, relates the transport properties of gauge 
theory plasmas to black hole dynamics. There are two 
transport coefficients of relativistic uncharged hydrody- 
namics at the first viscous order: the shear viscosity i] 
and the bulk viscosity C,. If there are charges present 
in the system, there are additional conductivities at this 
order. The gravitational calculation of the ratio of the 
shear viscosity to the entropy density s famously yields 
f = 3^ 01 ■ Low values of this ratio are rather generic in 
strongly coupled gauge theories, and seem to characterize 
the QCD plasma in heavy ion collisions. 

In Q the black hole dynamics was used to derive a 
simple formula for the ratio of the bulk viscosity to the 
shear viscosity: 






hH 



kH 



ds 



dp"- 



(1) 



where p" arc the charge densities associated with the bulk 
gauge fields, and (f)^ are the values of the bulk scalar 
fields on the horizon. The derivation employs the null 
horizon focusing equation, and uses the equilibrium black 
hole solution. When calculating the derivatives in ([T]), 
one keeps all the parameters of the theory such as cou- 
plings and masses fixed. This, as demonstrated in [3,Q, 
may require in general the knowledge of the complete 
thermal equilibrium gravitational background, and not 
only the horizon data. 

A natural question that arises is how transport coeffi- 
cients, and in particular the bulk viscosity, depend on the 
holographic radial direction r of the gravitational back- 
ground, which in the field theory translates to a depen- 
dence on the renormalization group (RG) scale. Trans- 
port coefficients depend on the couplings of the underly- 
ing microscopic field theory. Since the couplings in non- 
conformal field theories depend generically on the RG 
scale, it follows that the transport coefficients also have 



this property. Does this mean that the formula ([l} has an 
approximate regime of validity, for instance only at high 
temperature? There is evidence suggesting that this is 
not the case, and in fact it is the exact formula for the 
bulk viscosity Q-Ql . Still, a formal proof of its exactness 
is lacking. 

In this letter we will prove a much stronger statement 
in the framework of the fluid/gravity correspondence. We 
will study the dynamics of holographic screens (hypersur- 
faces) located in the radial direction at general r = r^, 
where tq is a constant. We will prove that all these 
holographic screens evolve, to all orders in the derivative 
expansion including all higher curvature corrections, ac- 
cording to the same hydrodynamics equations with iden- 
tical transport coefficients. The analysis will be carried 
out for normal fluids as well as for superfluids. An imme- 
diate consequence of this result is a proof of the exactness 
of the bulk viscosity formula ([T]) . 

Wc will consider [d + l)-dimensional gravitational 
backgrounds holographically describing thermal states in 
strongly coupled d-dimensional field theories. The (d+l)- 
dimensional gravitational action is / = Igravity [§] + 
l7natter[g,(l>i,A], where 
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where TZ is the Ricci scalar, A a (negative) cosmological 
constant and higher curvature refers to general higher 
curvature terms. We will set h = c = Gj^^ ~ 1 through- 
out. 

The matter action /,„a(ter[3, 0i,^] depends on (possi- 
bly charged) scalar fields 0^ and an (abelian) gauge field 
Ab- Wc assume that it has the general form 
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where Fab = 25[^Ab], Db = Vs — ieAB and Db = 
V B + ieAB for a given charge e. The Vi represent 




the potentials for the scalar fields. The field equations 
of ^ have solutions that holographically describe non- 
conformal normal fluids as well as superfluids. 

In this paper we will make use of two important con- 
servation laws. The first is that 
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This follows from the diffeomorphism invariance of the 
total action / and imposing the field equations of the 
matter (scalars and gauge field). The second is that 
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which follows from the gauge invariance of the matter 
action and imposing the scalar field equation. 

On-shcU, the metric field equations have the form 
Eab = 0, or more specifically, 



Eab — Gab + E[ab — T' 



AB 
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Gab and T™J** are the Einstein and matter stress-energy 
tensors respectively, while Hab represents any higher 
curvature corrections to the field equations. The field 
equations that arise from the variation of the matter ac- 
tion ^ are 
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where J is a source current that exists when e ^ 
and arises from the coupling between the scalars and the 
gauge field. The field equations of the scalars will be 
unimportant to our discussion. 

We will first consider solutions to this system that 
are dual to normal fluids on the boundary. In this 
case we take the scalar fields to be neutral e = and 
real, normalize the scalar and gauge field kinetic terms 
by setting Vi = l,y2 = | and have V3 an arbitrary 
potential for the scalar fields [3l- In (radially shifted) 
Eddington-Finkelstein coordinates the metric solutions 
can be parametrized with the following metric ansatz 



gABdx dx = dS(o) = -/( 



r + R)u^u,jdx^dx'^ 
2q{r + R)u^dx^'dr + h{r + R)P^^dx^'dx'' 
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Here we denote the coordinates of the bulk spacetime 
by a;^ = {r,x^),A = 0,...,d. The x^ can be thought of 
as local coordinates on hypersurfaces of constant r. At 
r = 0, the function / vanishes and there is a horizon. The 
vector u^ is a boost vector and has the dual role as a fluid 



velocity and the null normal to the horizon at r = 0. The 
indices /i, v etc. are raised with the flat metric 77^,^, and 
P^iu ~ t]fj,u + Ufj^Ui, is the projection tensor. The entropy 
density associated with this solution is s = jh{R)^^ . 

The function q can be chosen to be a constant by a 
gauge choice corresponding to a redefinition of the radial 
coordinate r. Thus, we expect that it will not play a 
role in the holographic description of the fluid. Another 
way to sec it, is to count the number of fields in the 
hydrodynamic description. In the absence of charges we 
have the energy density e, the pressure p and the fluid 
velocity u^, which is the same number of actual fields in 
(fTO| . namely /, h and w^. 

In addition to the metric, there are the solutions to the 
scalar fields (pi (r -I- R) and the gauge field 



Asdx^ = A{r + R)u^dx^ . 
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The one function in the gauge field is in line with the one 
extra charge p needed in the hydrodynamic description 
of a charged fluid. 

Let us allow now the variables, along with R and w^ 
to be functions of x'^. The zeroth order metric ([TO)) and 
other flelds as functions of x^ are no longer solutions of 
the field equations and must be corrected order by order 
in a derivative expansion. The fiuid/gravity correspon- 
dence proposes that at the nth order in the derivative ex- 
pansion [8|, the bulk constraint field equations projected 
on the time- like boundary at spatial infinity (r —^ 00), 
whose unit space- like normal is TVs, 
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are the relativistic Navier-Stokes equations c^'^T^""^ = 
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is the boundary stress- 
Similarly, in the presence of gauge fields. 



at the nth viscous order 

energy tensor 

the Gauss law constraint projected on the time-like 

boundary at infinity 
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is the additional current conservation law dfj,Jf^^^^ = in 
charged hydrodynamics, with J^^^-^^ the boundary current 
[9|, llO| . Alternatively, one can obtain these hydrodynam- 
ics equations from the Gauss-Codazzi equations of the 
horizon [ll|,|Il. 

Consider now how the constraint equations depend on 
r, that is what happens when we project them on an ar- 
bitrary time-like hypersurface r = tq. Suppose we have 
a gravitational solution at the (n — l)th order in deriva- 
tives, where n > 1, that is the bulk field equations are 
satisfied up to the (n — l)th order. We would like to 
analyze the structure of Eab at nth order. To impose 
the conservation law ([5]) at nth order we require that all 
the matter equations at this order satisfied. Then the 
conservation law takes the form 



^(^O)^AB — ^ ' 
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where the covariant derivative is with respect to the back- 
ground zeroth order metric (|10p and involves only radial 
derivatives. The fi components of this equation read 
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Using the unit normal to the hypersurface r = rg, iV — 
-^, N'^ = — i=u'^, we can recast ^ as 



9 f^Tp{n).TA^ 9 



(16) 

Solving, we find that (off-shell) the constraints have the 
form 






(17) 
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where i^^ (x^) is some arbitrary function of x''. 

Taking the limit ro — >■ oo, and renormalizing by an 
overall function of r, we get that Fi")(x'") - ^^'T|};^''^ = 
on-shell at the nth order. The factorization ([T7|) and the 
fact that the metric functions / and h are non-vanishing 
between the horizon and infinity lead to an important 
result: the bulk constraint equations projected on any hy- 
persurface r — To are the same hydrodynamics equations 
with identical transport coefficients ^. 

In the limit tq — > 0, \/f vanishes. At lowest orders in 
the derivative expansion, this is the near- horizon limit. 
It is well-known that in the limit as a "stretched horizon" 
r = To approaches the true horizon, ^N becomes the 
horizon normal u^^ 1141. Thus, we have 
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In the case of Einstein gravity, the right hand side is 
the Gauss-Codazzi equations for the null hypersurface, 
one of which is the null focusing equation. Therefore, 
following the derivation of ([1]) in [J|, we see that the ra- 
tio of the bulk viscosity to shear viscosity calculated at 
the horizon must be identical to the ratio that will be 
obtained from the boundary stress-energy tensor. Re- 
markably, the horizon calculation of bulk viscosity from 
the focusing equation requires only the zeroth order so- 
lution, while for the computation of the boundary stress, 
one needs to have the bulk solution to first viscous or- 
der. The bulk viscosity is an important transport co- 
efficient of field theories. The exactness of the formula 



^ This general result has been verified explicitly in conformal hy- 
drodynamics up to the first viscous order in \l3i . 



([T]) makes it a valuable tool in the study of strongly cou- 
pled non-conformal gauge theories and phenomenological 
holographic models of QCD. 

A similar result for the current conservation equation 
in charged hydrodynamics can be obtained using the 
Gauss law constraint. As before, assume that we have 
a solution to the gauge field equations (n— l)th order, so 
that the first non-zero piece of E^ is at nth order. We 
must also impose the scalar field equations at nth order. 
Then the conservation law ([71) reduces to 



a^(g/iV£;{»)'-) = 



The constraint (fT3l) at nth order has the form 
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Thus, we have an equation for the constraint, which 
yields 



h'^Vf 



(21) 



for some arbitrary function G^"' (x^), which on-shell must 
be proportional to the boundary current conservation 
equations. The same analysis as above applies and this 
constraint is the same current conservation equation on 
each r = rQ slice including the horizon, with identical 
transport coefficients. 

Consider now solutions to the general gravitational sys- 
tem that holographically describe relativistic supcrfiuids. 
In this case, the solutions can be parametrized with the 
following metric ansatz [l5l - ll8l | 

(is?Q-) = —f{r + R)u^Ui,dx^dx'^ 

+ 2q{r + R)Ufj^dx^dr + k{r + R)n^n^dx''dx'' 

-\- j{r -\- R){Ufini, + Uvnij)dx^dx'^ + r^P^^dx^dx'^ . 

(22) 



Here n^ satisfies u^n,, 



0, n^n. 



1 and P, 
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The ansatz for the gauge field solution 



takes the form 



Asdx^ = A{r + R)Ufj,dx'' 



B(r 



R)nfj_dx^ 
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The five functions in the metric and gauge field, together 



with u^ and 



match the number of fields needed to 



characterize the relativistic superfiuid. In addition to 
the fields associated with the normal fluid (e,p, p, u**), we 
have a new vector field n^, which is the component of the 
superfiuid velocity orthogonal to the normal fiuid velocity 
u^, the superfiuid energy density eg and the superfiuid 
charge density ps- Thus, as before, we expect that we 
can gauge away the function q. 

We can now repeat the analysis of how the conservation 
laws effect the constraints for this ansatz. The result 



is a similar factorization of the constraints into a radial 
function times a function of x'^ , 
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(24) 
(25) 



Therefore, we conclude that in the superfluid case there 
is again no running of the transport coefficients once the 
UV boundary conditions have been fixed. 

It is straightforward to check that the inclusion of a 
Chern-Simons term in the bulk action, corresponding to 
having quantum anomalies in the dual field theory, docs 
not change the above discussion. Thus, the parity violat- 
ing anomalous transport coefficients in the hydrodynam- 
ics are independent of the radial direction, as expected 
from the field theory analysis for normal fluids |19l423 | 
and superfluids [11 iM Hi ■ 

The above analysis applies also to the hydrodynamic 
expansion around Rindlcr space [ll|, |25| , where the trans- 



port coefficients depend on the UV boundary conditions 
at the cutoff hypersurface r = Tc, but do not depend ex- 
plicitly on the holographic radial direction r from the UV 
to the IR [13. 

Finally, note that the study in this letter is different 
than the type of KG flow discussed in [23| . In these works 
one imposes boundary conditions at a cutoff hypersurface 
and finds the corresponding (non-singular) interior bulk 
solution. Changing r = r,. changes the solution and its 
corresponding energy-momentum tensor, leading to a RG 
flow in the transport coefficients. In our case, we consider 
the hypcrsurfaccs rg in given solution. 
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